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We study the Quantum Electrodynamics of 2D and 3D Dirac semimetals by means of a self- 
consistent resolution of the Schwinger-Dyson equations, aiming to obtain the respective phase di¬ 
agrams in terms of the relative strength of the Coulomb interaction and the number N of Dirac 
fermions. In this framework, 2D Dirac semimetals have just a strong-coupling instability character¬ 
ized by exciton condensation (and dynamical generation of mass) that we find at a critical coupling 
well above previous theoretical estimates, thus explaining the absence of that instability in free¬ 
standing graphene samples. On the other hand, we show that 3D Dirac semimetals have a richer 
phase diagram, with a strong-coupling instability leading to dynamical mass generation up to N = 

4 and a line of critical points for larger values of N characterized by the vanishing of the electron 
quasiparticle weight in the low-energy limit. Such a critical behavior signals the transition to a 
strongly correlated liquid, characterized by noninteger scaling dimensions that imply the absence of 
a pole in the electron propagator and are the signature of non-Fermi liquid behavior with no stable 
electron quasiparticles. 


I. INTRODUCTION 

The discovery of graphenejlj has marked the beginning of a new chapter in condensed matter physics, with the 
appearance of new fundamental concepts and materials with unconventional properties. We have known about other 
genuine 2D materials like the transition metal dichalcogenides, and we have learned about the features hidden in the 
band structure of the topological insulators (3,0- Lately, wc have also seen the discovery of 3D Dirac semimetals which 
are the higher-dimensional analog of graphene U{Z], and we have witnessed the ongoing search of Weyl semimetals 
with a built-in breakdown of parity and time-reversal invariance. 

In the above instances, most part of the unconventional features of the materials come from the peculiar geometrical 
and topological properties of the band structure. Moreover, these are electron systems that are prone to being placed 
in the strong coupling regime, with a large relative strength of the Coulomb interaction. Graphene should be a clear 
example of electron system with strong interaction, given the large ratio between the square of the electron charge 
and the Fermi velocity in the 2D material. However, graphene is not a prototype of strongly correlated system, even 
in the case of the free-standing material in vacuum, as experimental observations have shown no sign of electronic 
instability down to very low doping levels 

On the other hand, most part of theoretical studies [9l42l| have estimated that 2D Dirac semimetals should have an 
excitonic instability at a critical point below the maximum interaction strength attained in graphene suspended in 
vacuum. The absence of any signature of a gap in the electronic spectrum, even below the meV scale, is certainly a 
puzzling evidence regarding the behavior of the 2D material. Given that the theoretical analyses have been mainly 
based on a ladder approximation to the electron self-energy corrections, the discrepancy between theory and experi¬ 
ment calls into question the use of such approximate methods in electron systems that are placed in the strong-coupling 
regime [22j. 

In this paper, we apply a nonperturbative approach to the investigation of the effects of the Coulomb interaction in 
both 2D and 3D semimetals, with the aim of mapping more confidently the different phases that may appear in those 
electron systems. More precisely, we carry out the self-consistent resolution of the Schwinger-Dyson equations for the 
Quantum Electrodynamics (QED) of 2D and 3D Dirac semimetals, in which the scalar part of the electromagnetic 
potential is used to mediate in each case the long-range e-e interaction. This approach has to be implemented in general 
with some kind of truncation to guarantee its practical feasibility. In this regard, we have relied on a formulation 
of the equations that amounts to including all kinds of diagrammatic contributions except those containing vertex 
corrections. Nevertheless, the solutions obtained in this way account for the renormalization of all the quasiparticle 
parameters, giving rise to frequency and momentum-dependent forms of the quasiparticle weight, the Fermi velocity 
and the dynamical Dirac fermion mass. 

In this framework, we will see that 2D Dirac semimetals have just a strong-coupling instability characterized by 
exciton condensation (and dynamical generation of mass) that we find at a critical coupling well above the estimates 
based on a ladder approximation, thus explaining the absence of that instability in free-standing graphene samples. 
On the other hand, we will show that 3D Dirac semimetals have a richer phase diagram, with a strong-coupling 
instability leading to dynamical mass generation up to TV = 4 and a line of critical points for larger values of TV 
characterized by the vanishing of the electron quasiparticle weight in the low-energy limit [23j] . We will see that such a 
critical behavior marks the transition to a strongly correlated liquid, characterized by noninteger scaling dimensions 
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that imply the absence of a pole in the electron propagator and are the signature of non-Fermi liquid behavior with 
no stable electron quasiparticlesQ- 


II. QUANTUM ELECTRODYNAMICS OF DIRAC SEMIMETALS 


We focus on the QED of Dirac semimetals, for which the Fermi velocity Vp is much smaller than the speed of light. 
The dynamics of these systems can be then described by the interaction of a number N of four-component Dirac 
spinor fields ^(r) representing the electron quasiparticles and the scalar part c/>( r) of the electromagnetic potential. 
In principle, each spinor can represent the electronic states around a different Dirac point in momentum space, but 
we will not make more explicit such a discrimination as it does not play any role in the subsequent analysis. The 
hamiltonian can be written in general as 


H = ivp 


d D r ipj (r) 7 o 7 • Vipi( r) + e 


d D ripl(r)ipi(r) </>(r) 


( 1 ) 


where { 7 qJ is a set of Dirac matrices satisfying { 70 , 7 / 3 } = 2 ? 7 Q /3 (with the Minkowski metric in D + 1 dimensions 
77 = diag(-l,l,...l)). 

The expression m holds equally well for dimension D = 2 and 3, but the propagator of the </> field is very different 
in the two cases. The scalar field has to mediate the e-e interaction with long-range Coulomb potential F(r) = l/|r|, 
irrespective of the spatial dimension. At D = 2, this leads to a bare propagator D 0 (q, oj) for the (j) field 

A)(q)| D=2 = (2) 

while in 3D space the bare propagator is instead 

-Do(q)|_o = 3 = ^2 (3) 

The effects of the interaction can be characterized through the corrections undergone by the scalar and the Dirac 
field propagators. The full Dirac propagator G(k, ujk) has in general a representation of the form 

G(k, w fe ) -1 = (w fc - uf7o 7 ' k ) - s (k. Wfc) (4) 


in terms of a self-energy correction E(k, Wfc) that contributes to renormalize the bare quasiparticle parameters. This 
object is given in turn by the equation 


*E(k, w fc ) 


d D p doj p 
(2 n) D 27r 


D( p, w p )G(k - p, a;*, - w p )r(p, k, w fc ) 


(5) 


where D(p,ui p ) stands for the full propagator of the scalar potential and T(q, k, w*,) represents the irreducible 
three-point vertex. More precisely, this function is defined by the expectation value 

zer(q,w g ;k,w fc ) = (</>(q, w g )V' l (k - q,w fc - w fc ))i P i ( 6 ) 


where 1PI means that we must take the irreducible part of the correlator. 

Furthermore, D( p,w p ) has also its own equation representing it in terms of the full propagators and the irreducible 
vertex. We can write 


D( q,w g ) 1 = D 0 (q) 1 -n(q,w g ) 


(7) 


with the polarization II(q, ui q ) being given by 


*n(q, LOq) = Ne 2 


d D p dujp 

(2n) D ^~ 


Tr [G(q + p, u q + w p )r(q, p, w p )G(p, w p )] 


( 8 ) 


Moreover, the form of the irreducible vertex r(q, u> q \ k, u>k) is constrained by the Ward identity aris ing from the 
reduced gauge invariance of the model, admitting also a representation in terms of the full propagators |25j . 

The expressions (JoJ) and © correspond to the Schwinger-Dyson equations of the model. Together with a suitable 
representation of the irreducible vertex, they may lead to valuable information about the form of the full fermion and 
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interaction propagators. In general, however, one has to resort to some kind of truncation to achieve a self-consistent 
resolution of the integral equations. In what follows, we will apply a common procedure, the so-called bare vertex 
approximation, by which r(q, u q \ k, u>k) is set equal to the unit matrix in the resolution of ([5]) and ©. We note that 
this truncation does not satisfy the mentioned Ward identity, which relates the irreducible vertex to the derivative 
with respect to the frequency of the fermion self-energy. In this regard, the present work focuses on the investigation 
of dynamical effects in the renormalization of the quasiparticle parameters, which could be further improved by 
introducing a suitable ansatz for the vertex (in similar fashion as in the fully relativistic QED^). Nevertheless, we 
note that the relevant features re por ted below within our approach are consistent with the results found by means 
of renormalization group methods [23|, which supports the present formulation. The main advantage of adopting the 
mentioned truncation is that it leads to a very convenient implementation of the self-consistent approach, allowing us 
to attain easily convergence in the recursive resolution of the Schwinger-Dyson equations. 

Without the vertex corrections, d5j) and © lead indeed to closed self-consistent equations, shown diagrammatically 
in Fig. □ This representation allows us to establish a comparison with other standard approaches used to deal with 
many-body corrections. In particular, it becomes clear that the contributions accounted for by the diagrams in Fig. 
□ have a much more comprehensive content than other approaches dealing with the RPA sum of bubble diagrams for 
the polarization. This makes the present computational scheme much more reliable to describe the electron system 
away from the weak-coupling regime, incorporating effects like the renormalization of the Fermi velocity and the 
quasiparticle weight which are essential to capture the different critical points of the Dirac semimetals. 




FIG. 1: Diagrammatic representation of the Schwinger-Dyson equations (in the bare vertex approximation). The 
thick(thin) straight line represents the dressed(free) Dirac fermion propagator and the thick(thin) wiggly line represents the 
dressed (undressed) interaction propagator. 


III. SELF-CONSISTENT RESOLUTION OF SCHWINGER-DYSON EQUATIONS 

The integral equations represented in Fig. |T] can be solved numerically by applying a recursive procedure, after 
rotating first all the frequencies in the complex plane, u; = — ico, to make the passage to a Euclidean space in 
the variables (uJ, k). In practice, the integrals can be done numerically by discretizing the frequency and momentum 
variables. By choosing a set of frequencies oJ = 7r(2n+l)T with n = 0, ±1, ±2,..., we can interpret such a discretization 
as the result of placing the theory at finite temperature T. On the other hand, computing with a grid in momentum 
space is equivalent to describing a system with finite spatial size. In this case, we can check the finite-size scaling of 
the results in order to extrapolate the behavior over large distances. 

For the self-consistent resolution, it becomes convenient to represent the fermion propagator in terms of renormal¬ 
ization factors z^( k, iUJ) for the electron wave-function and z v ( k, iTU) for the Fermi velocity, adding moreover another 
factor z m (k.,iU) to allow for the dynamical generation of a mass for the Dirac fermions. Thus we write the full Dirac 
propagator in the form 

G(k, iui) = (z^( k, iuj)iuJ - z v ( k, k - z m ( k, iUJ) 7o ) _1 ( 9 ) 

In this way, the resolution consists in finding the functions 2^(k, iUJ), z„(k, ioj) and z m (k,iuJ) that attain the self- 
consistency in the Schwinger-Dyson equations. 

One more important detail is that the polarization II(q, uj q ) may develop spurious divergences when computing the 
momentum integrals with a simple cutoff A^. In general, only a g aug e-invariant regularization scheme can produce 
results without non-physical power-law dependences on the cutoff[26j. These are anyhow additive contributions to 
the polarization, which makes possible to get rid of them by a suitable subtraction procedure. Thus, computing with 
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the momentum cutoff, the polarization at D = 2 shows a contribution proportional to A*,, while the corresponding 
function at D = 3 has terms growing as large as A 2 . In our self-consistent resolution, we have carried out the 
frequency integrals first with a cutoff A upA/., implementing afterwards the subtraction procedure to remove the 
power-law dependences on A*, from the polarization. In this way, we have ended up with expressions of II(q, uj q ) that 
are functionals of the renormalization factors, displaying leading behaviors at small q proportional to |q| and q 2 , 
respectively, for D = 2 and D = 3. 


A. 2D Dirac semimetals 


Solving the Schwinger-Dyson equations at D = 2, we find in general that the function z^( k, iuj) giving the quasipar¬ 
ticle weight remains bounded, while z v ( k, iuj) diverges in the limit of small momentum k — > 0. As long as the effective 
Fermi velocity depends on the momentum scale, it is convenient to define the bare value vb = z v ( A^, 0)up, which can 
be taken as a good measure of the Fermi velocity at the microscopic scale (it is always verified that z^{ Afc,0) ss 1). 
We can then define the unrenormalized coupling giving the bare interaction strength as 

a = e 2 /4:irvB ( 10 ) 


which can take different values depending on the particular Fermi velocities of the 2D Dirac semimetals. 

As an illustration of the general behavior, Fig. [^represents the solution obtained for z^{ k, iuj) and z v ( k, iuj) for N = 
2 and a = 2.2, that is, for parameters that should be appropriate to describe graphene samples suspended in vacuum. 
The resolution has been carried out taking a discretization of the frequency variables such that 2nT ps 0.01 eV. In 
this case, the self-consistency in the equations is only attained when z m (k, iuj) is set identically equal to zero. The 
behavior found for z^,(k, iuj) and z v ( k, iuj) is in agreement with the general trend obtained from renormalizatio n group 
methods, which found the divergence of the Fermi velocity in the low-energy limit as a most relevant feature[26l 1271]. 

For comparison, we have also represented in Fig. [2] the dependence on the energy scale e of the inverse of the 
quasiparticle weight z(e) and the renormalized Fermi velocity v(e) obtained from the renormalization group approach 
in the large-./V approximation[27|• It can be observed anyhow that the plot of z v (k, 0)/zj,(k, 0) (full line in Fig. [21(b)) 
follows the experimental results of Ref. HI (Fig. 2(c) in that paper) much more accurately than the scale dependence 
of the Fermi velocity obtained with the renormalization group method in the large-iV approximation (dashed line in 
Fig. Hb))[13. 
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(a) 


(b) 


FIG. 2: Plot of the factors z^,(0,iu}) and z v (k, 0)/x;^(k, 0) (full lines in (a) and (b)) for a 2D Dirac semimetal with N = 2, bare 
coupling a = 2.2, and 2ttT ~ 0.01 eV. The dashed lines represent the dependence on the energy scale e of the inverse of the 
quasiparticle weight z(e) (in (a)) and the renormalized Fermi velocity u(e) (in (b)) obtained with the renormalization group 
approach for the same bare coupling in the large-./V approximation. 

With our nonperturbative approach, moreover, we can ask whether the tendency towards a noninteracting Fermi 
liquid, implied by the growth of the Fermi velocity, can be arrested by some instability as the bare coupling a is 
increased. The outcome of this search is that the other relevant feature of the 2D Dirac semimetals is the development 
of a nonvanishing mass z m (k,iuj) at sufficiently large interaction strength, as illustrated in Fig. OE This corresponds 
to the onset of a phase with chiral symmetry breaking and dynamical generation of a gap for the Dirac fermions, as 
predicted by several other methods [9l-[2l|. 

We observe that the divergence of z v in the low-energy limit does not prevent the development of a nonvanishing 
dynamical mass z m , while z,p remains finite accross the transition. This latter fact implies that chiral symmetry 
breaking proceeds without an anomalous scaling of the Dirac fermion field, in agreement with field theory studies of 
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FIG. 3: Plot of 2^,(0, iio), 2 „(k, 0)/2,/,(k, 0), and 2 m (k, 0)/2,/,(k, 0) measured in eV (full lines in (a), (b) and (c)) for a 2D Dirac 
semimetal with N = 2, bare coupling a = 3.37, and 2nT ~ 0.01 eV. The dashed lines represent the dependence on the energy 
scale e of the inverse of the quasiparticle weight z(e) (in (a)) and the renormalized Fermi velocity v(e ) (in (b)) obtained with 
the renormalization group approach for the same bare coupling in the large-IV approximation. 


that phenomenon H^]. In general, we also expect that the divergent growth of the Fermi velocity is arrested at the 
energy scale for which the quasiparticle dispersion becomes gapped (which is not appreciated in the case of Fig. |3](b) 
as the infrared cutoff set in that plot by finite-size effects is slightly above 1 meV). 

The present approach has the virtue of allowing an accurate determination of the critical coupling a c for the 
transition to the broken symmetry phase. We have represented in Fig. [J]the plot of the critical coupling obtained 
as a function of N , which leads to a map of the two different phases in the QED of the 2D Dirac semimetals. In 
agreement with earlier analyses, we observe that a c turns out to grow with N, though in the present resolution there 
seems to be no upper limit in the number of Dirac fermions for the development of the transition[30j|. We have also 
checked that the approach to the critical coupling seems to be consistent with a transition of infinite order, since the 
dynamical mass exhibits an inflection point as a function of coupling constant above a c which is the signature of that 
kind of transition under finite-size effects [31] . 



FIG. 4: Phase diagram of the QED of 2D Dirac semimetals showing the region with dynamically generated fermion mass Zm ^ 0. 
The critical line has been obtained from the self-consistent resolution of the Schwinger-Dyson equations with 2nT « 0.01 eV. 
The dashed line marks the nominal interaction strength a « 2.2 corresponding to graphene suspended in vacuum. 

In any case, the most relevant result regarding the phase diagram in Fig. [I] is that the point corresponding to real 
graphene samples suspended in vacuum falls in the region with no dynamical generation of mass. For the number 
IV = 2 of Dirac fermions in graphene, the critical coupling obtained for chiral symmetry breaking turns out to be 
indeed well above most part of previous estimates relying on a restricted sum of many-body corrections. The present 
results explain therefore that no gap has been found in the electronic spectrum of graphene, even in experiments 
looking very close to the Dirac point [§]. The reason for the unexpectedly large values of the critical coupling can be 
traced back to the combination of the slight suppression of the quasiparticle weight and the large growth of the Fermi 
velocity at low energies [H]. These two effects cooperate to reduce significantly the effective strength of the Coulomb 
interaction for the development of the gap, stressing the importance of a proper account of all the renormalization 
factors for the accurate determination of the transition to the broken symmetry phase. 
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B. 3D Dirac semimetals 

The 3D Dirac semimetals have in general a number of Dirac points that have attached (each of them) fermions with 
the two different chiralities. This is in particular the case of materials recently discovered like NaaBi and CdaAs 2 , as 
already clarified by their theoretical analysis [HI [34|. Such a distinctive feature of the 3D Dirac semimetals is relevant 
in the present study, since it makes possible the dynamical generation of mass and opening of a gap in these systems 
from the hybridization of two chiralities at the same Dirac point. 

More precisely, the low-energy electronic states in both NasBi and CdaAs 2 can be naturally arranged into four- 
component spinors around each of two Dirac points, in such a way that the Dirac matrices appear in the chiral 
representation 


7o7i = m ( q 1 -<j t ) ( n ) 

with 771 = 772 = 1 and 773 = ±1 depending on the Dirac point. The dynamical mass generation that we impose with 
the ansatz ([9j corresponds in this scheme to the mixing of the two chiralities, realized by the Dirac matrix 

70 = ( \ - 0 * ) (12) 

A term proportional to (fl2l) has been identified in Refs. (33) and [34j as one of the possible perturbations of the 
Dirac hamiltonian in NaaBi and Cd 3 As 2 - The physical meaning of such a term has to be found in the breakdown of 
the threefold and fourfold rotational symmetry in each case, having an effect similar to that induced by strain in the 
crystal lattice. 

The dynamical mass generation we are discussing corresponds then to the development of an expectation value 

($• 7o^) t^O (13) 

within each Dirac point i. We can imagine nevertheless the possibility of a symmetry breaking pattern with order 
parameter given by 


{ipjMipj) ± 0 (14) 

with a suitable matrix M and a pair of Dirac points i ^ j [§5[. Such a condensation has the feature of involving a 
finite momentum Q, needed to connect different Dirac points. This can be pictured diagrammatically, as the order 
parameters m and (gH) can be characterized in terms of the respective three-point vertices in Fig. [5j The point 
is that, as implied by the renormalization group approach in Ref. [36( for the 2D case, one has to hybridize the two 
chiralities at the same Dirac point so that the long-range Coulomb interaction may induce the strongest scaling of 
the three-point vertex, leading eventually to the condensation signaled by m■ The vertex in Fig. [5jb) is built from 
spinors at different Dirac points, which will not map in general onto each other upon a rigid shift by the momentum 
Q and may lead therefore to a weaker overlap. This means that, for a generic 3D Dirac semimetal with long-range 
Coulomb interaction, the maximum strength will be set by the vertex in Fig. [U(a), ensuring at least the condensation 
given by (1131) whenever symmetry breaking is to take place in the system. 



(a) 


(b) 


FIG. 5: Diagrammatic representation of the three-point vertices involving the composite operators i/dyoi/ji and ipjMilij, with 
respective momentum transfer q = 0 and q = Q connecting Dirac points. 

Apart from the phenomenon of dynamical mass generation, the 3D Dirac semimetals have a tendency to develop at 
strong coupling a drastic attenuation of the quasiparticle weight at low energies, with a much softer renormalization 
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of the Fermi velocity in comparison to their 2D Dirac analogues (37}. This behavior has been already found in an 
analytic study of the 3D electron systems in the large -N limit, where it has been possible to establish rigorously the 
existence of a critical coupling at which the quasiparticle weight vanishes in the low-energy limit[23|. In Fig. [ 6 ] we 
can see the effect of the critical behavior in the functions 2 ^,(k, iUJ) and z v (k,iUJ), computed in the present approach 
for IV = 6 close to the critical point and with a discretization such that 2ttT ss 0.02 eV. For smaller values of N, 
we will see that there is however an interplay between that quasiparticle attenuation and the tendency to dynamical 
generation of mass. For N < 4, this latter effect becomes actually dominant, leading to a phase with chiral symmetry 
breaking that is the analog of the broken symmetry phase found in the 2D Dirac semimetals |38j|. 



FIG. 6 : Plot of the factors z^(0,iu}) and z„(k, 0)/z q p(U., 0) for a 3D Dirac semimetal with N = 6, bare coupling g = 36.8, and 
2 t tT « 0.02 eV. 


From the self-consistent resolution of the Schwinger-Dyson equations, we have determined for each value of N the 
critical coupling at which the electron system becomes first unstable, either from the vanishing of the quasiparticle 
weight or from the dynamical generation of a gap. In the case of the 3D Dirac semimetals, one has to take special 
care to refer the parameters to a given scale, since quantities like the electron charge and the Fermi velocity are 
renormalized at low energies. In this respect, we have chosen to define the bare electron charge es at the highest 
value of the momentum cutoff, according to the relation e 2 B = A^Z?(Afe, 0)e 2 . Then, we can take for the microscopic 
parameter e B the standard value of the electron charge. As in the case of the 2D Dirac semimetals, we have also 
defined the bare Fermi velocity by v B = z v (Ak,0)vF ■ Thus, we have computed all the critical couplings referred to 
the relative interaction strength at the microscopic scale, given by 

g = Ne 2 B /4TTV B (15) 



FIG. 7: Phase diagram of the QED of 3D Dirac semimetals, showing the region with dynamical generation of mass ( z m 7 ^ 0) 
and the region corresponding to non-Fermi liquid behavior (with Re(^(0, iui)) oc a; 7 , 7 < 0). The critical lines have been 
obtained from the self-consistent resolution of the Schwinger-Dyson equations with 2nT tts 0.02 eV. 

The results we have obtained are condensed in the phase diagram shown in Fig. [7] We observe that there is always 
a phase connected to weak coupling for all values of N, characterized by a gapless spectrum of quasiparticles whose 
parameters remain regular at low energies. This phase terminates for N < 4 in the dynamical generation of a fermion 
mass at sufficiently strong coupling, which in our approach is reflected in the onset of a nonvanishing z rn (k, ioj). This 
is shown in Fig. [ 8 l where we have represented the different renormalization factors for N = 2 at a coupling above 
the critical point. For this value of N, we observe that the renormalization of the quasiparticle weight is a moderate 







effect as the gap opens up in the electronic spectrum. This soft behavior has been also observed in the studies of 
dynamical mass generation in the fully relativistic QED, carried out by means of the self-consistent resolution of the 
corresponding Schwinger-Dyson equations [25|. 





a> k k 


(b) 


(c) 


FIG. 8: Plot of the factors z^,(0, iu>), z v (k, 0)/^(k, 0) and z m (k, 0)/s^(k, 0) (measured in eV) for a 3D Dirac semimetal with 
N = 2, bare coupling g = 15.1, and 2nT « 0.02 eV. 


The case of N = 4 is however specially interesting, since there is then an interplay between the dynamical generation 
of mass and the strong attenuation of the quasiparticle weight. This can be observed in Fig. [9j where we have plotted 
ico), z v (k,iui) and 2 m (k, iUJ) for different couplings below and above the point where the mass develops. We see 
that, while the breakdown of chiral symmetry takes place before the system is completely destabilized by the large 
growth of 2^(k, iw), this latter effect may still have a large impact on the observation of the quasiparticles in the 
electron system. 





lj k k 


(a) 


(b) 


(c) 


FIG. 9: Plot of the factors «^,(0, iui), z„(k, 0)/z,/,(k, 0) and 2 m (k, 0)/z^(k, 0) (in eV) for a 3D Dirac semimetal with IV = 4, 
bare coupling g = 30.1, 27.6, 24.9 (from top to bottom in (a) and (c), from bottom to top in (b)), and 2nT « 0.01 eV. 

On the other hand, we find that for N > 4 there is always a critical coupling at which the divergence of z^,(0,iuj) 
takes place before the dynamical generation of mass, according to the trend illustrated in Fig. [G] The present 
approach allows us moreover to investigate the phase of the electron system above the critical point. The self- 
consistent resolution of the Schwinger-Dyson equations gives rise in that case to renormalization factors that get in 
general an imaginary part, as shown in Fig. 1101 This has to be interpreted as the signature of a nonperturbative 
instability of the electron quasiparticles since, in the conventional perturbative approach, the self-energy corrections 
obtained after Wick rotation uj = ioJ can only account for the renormalization of the real part of the quasiparticle 
parameters. In our approach, the divergence of the imaginary part of z^(0,ioJ) at uJ = 0 points actually to the 
development of a strongly correlated liquid, which is confirmed by the concomitant suppression of the quasiparticle 
weight in the low-energy limit, observed in the plots of Fig. [10] 

A detailed analysis of the renormalization factor z^(0,iZJ) in Fig. [TOl reveals indeed that the real part of such a 
function follows accurately a power-law behavior as w -> 0. This can be clearly seen in the plots of Fig. |TT] where 
aJRe(^(0,ioJ)) is represented in linear and logarithmic scale. The fit to the power-law dependence 

ZJ Re(^(0, iuJ)) oc 


(16) 
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FIG. 10: Plot of the real and imaginary parts of the factors z^,(0,iu>) and z„(k,0) for a 3D Dirac semimetal with N = 6, bare 
coupling g = 37.2, 37.4, 37.6 (from lower to higher absolute value in both sides of the plot), and 2nT ss 0.02 eV. 


gives an exponent n ss 0.7, with little variation between the three curves for the different couplings. Equivalently, 
this corresponds to having a nonvanishing anomalous dimension of the Dirac fermion field with a value of ~ 0.3. The 
inspection of the renormalization of vf, dictated by the function Re(,z„(k, 0))/Re(^(k, 0)), leads in this case to a 
picture very similar to that in Fig. 161(b). This shows that /i corresponds to the exponent governing both the frequency 
and momentum scaling of the electron propagator (implying therefore a dynamical critical exponent equal to 1). 




CJ 0) 


(a) 


(b) 


FIG. 11: Plot of l o Re(z^,(0, iu))) for a 3D Dirac semimetal with the same parameters as in Fig. 1101 with (a) linear and (b) 
logarithmic scale (with the three different curves virtually collapsed onto the same line). 


From the physical point of view, the main consequence of the noninteger exponent /i is the absence of a pole in 
the electron propagator, which reflects the lack of low-energy fermion excitations. This is characteristic of correlated 
systems in low dimensions, where the interactions may drive into a non-Fermi liquid phase with nontrivial scaling 
exponents [39M43| . Our system provides in this respect an example of such a behavior at D = 3, illustrating moreover 
the transition to the strongly correlated phase from the renormalized Fermi liquid, which is the phase of 3D Dirac 
semimetals at sufficiently weak coupling. 


IV. CONCLUSION 

We have seen that the behavior of 2D and 3D Dirac semimetals is governed by quite different effects in their 
respective strong-coupling regimes. In our approach to the 2D semimetals, there is in general a tendency of the Fermi 
velocity of quasiparticles to grow large in the low-energy limit, in agreement with the renormalization group studies 
carried out in the large -N approximation(27j. In the case of the 3D Dirac semimetals, we observe instead that the 
electron system is prone to develop an attenuation of the quasiparticle weight, with a less significant renormalization 
of the Fermi velocity. This tendency has been also identified in the large -N limit of the 3D Dirac semimetals, which 
shows the existence of a critical coupling characterized by the divergence of the electron scaling dimension j23|. 

In the 2D Dirac semimetals, the divergence of the Fermi velocity in the low-energy limit is the dominant feature 
that explains for instance the absence of significant correlation effects in the graphene layer. Our nonperturbative 
solution of the Schwinger-Dyson equations incorporates naturally the scaling of the Fermi velocity, allowing us to 
reach very good agreement with the experimental measures from graphene samples at very low doping levels |8[. As a 
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result of such a renormalization, we have found that the interaction strength has to be set to relatively large values, 
above those attained in the suspended graphene samples, in order to open up a phase with exciton condensation and 
dynamical generation of a gap in the 2D Dirac semimetals. 

The picture changes into a richer phase diagram for the 3D Dirac semimetals, as a consequence of the interplay 
between the attenuation of the electron quasiparticles that prevails at large N and the tendency to dynamical mass 
generation (analogous to the chiral symmetry breaking of the fully relativistic OKI) Fill I77TI | 1 that is dominant at small 
N. Both effects seem to coexist at the interface found for a number of Dirac fermions N = 4. Most interestingly, 
our self-consistent resolution has also revealed the phase of the system above the large-TV critical point, allowing us 
to characterize the properties of a strongly correlated liquid that is reminiscent of other systems with suppression of 
electron quasiparticles making the transition from marginal Fermi liquid [39| to non-Fermi liquid behavior [40l443j | . 

We remark that our analysis of the 3D Dirac semimetals can be extended to map also the large -N regime of 3D Weyl 
semimetals. These are a class of semimetals in which a number of Weyl points host fermions with a given chirality, 
represented in terms of two-component spinors. This means that a self-consistent resolution of the Schwinger-Dyson 
equations may be also carried out for these systems, writing now the full propagator of the Weyl fermions around a 
given Weyl point as 

G(k, iZJ) = (^(k, — 2 t ,(k, iui)vFcr • k) 1 2 3 4 5 (17) 

We may parallel the above approach to predict the existence of a phase with suppression of electron quasiparticles, 
similar to that found for the 3D Dirac semimetals and covering the right part of the phase diagram in Fig. [7] At small 
N (taken now as the number of pairs of Weyl points), a strong-coupling phase corresponding to fermion condensation 
may also arise, with an order parameter mixing fermions at different Weyl points as in (1141) . The strength of this 
instability cannot be assessed generically, however, since it may be highly dependent on the particular overlapping 
of spinors from different Weyl points. Nevertheless, we may expect a strong-coupling symmetry broken phase for 3D 
Weyl semimetals at small N, with a phase boundary determined by the particular form of the spinors in the material 
hosting the Weyl points. 

We finally comment on the feasibility to observe the strong-coupling phases of the 3D semimetals, according to 
the values of the Fermi velocity and number of Dirac or Weyl points found in different materials. We recall in this 
regard that the best known examples of 3D Dirac semimetals (NasBi and CdaAs 2 ) have a number N = 2 of Dirac 
fermions and Fermi velocities that have been measured with certain accuracy pi Fhi . The quasiparticle dispersion 
shows in both cases an anisotropy that is reflected in the values of vf , which may be reduced by a factor of ~ 4 in 
one of the directions in momentum space [53j. We can make conservative estimates of the coupling defined in (fl5l) by 
taking the largest Fermi velocity in each case, with the result that g ~ 10 for NasBi and g ~ 3 for CdsAs 2 . These 
couplings turn out to be below the critical coupling for dynamical mass generation, which corresponds to the critical 
point at g* ~ 14.7 for N = 2 in the phase diagram of Fig. [7] This places the two mentioned 3D Dirac semimetals in 
the gapless phase, showing that the transition to the regime with dynamical mass generation would require a Fermi 
velocity at least about 50% smaller than the largest value in NaaBi. 

On the other hand, there should be good prospects to observe the strong-coupling phase at large N in Weyl 
semimetals. The most promising candidates for this class are the pyrochlore iridates and TaAs, which have 12 pairs 
of Weyl points. With the value of the Fermi velocity measured for TaAs[54|, the large value of N already sets the 
effective coupling g for this material well above the critical coupling for the non-Fermi liquid regime, which is at 
g* « 21.1 for N = 12. This leads us to expect that such a material should not behave as a regular Fermi liquid when 
observed at filling levels sufficiently close to the Weyl points. In general, we conclude that the strong-coupling phases 
that we have studied in this paper are not beyond reach, and that they may be found in systems with reasonably low 
Fermi velocities, or with a sufficiently large number of Dirac or Weyl points already exhibited by several materials. 

We acknowledge the financial support from MICINN (Spain) through grant FIS2011-23713 and from MINECO 
(Spain) through grant FIS2014-57432-P. 
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